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Abstract 



CN , Structure transformation by heat treatment in single-walled carbon nanotubes (SWCNT) is 

investigated using molecular dynamics simulation. The critical temperature for the collapse of 

"^ ' pure SWCNT is as high as 4655 K due to strong covalent carbon-carbon bonding. Above 2000 K, 

the cross section of SWCNT changes from circle to ellipse. The self-repairing capability is then 

investigated and two efficient processes are observed for the SWCNT to repair themselves. (1) In 

Ch ' the first mechanism, vacancy defects aggregate to form a bigger hole, and a bottleneck junction 

'^ ■ is constructed nearby. (2) In the second mechanism, a local curvature is generated around the 

O . isolate vacancy to smooth the SWCNT. Benefit from the powerful self-repairing capability, defective 

SWCNT can seek a stable configuration at high temperatures; thus the critical temperature for 

^ ■ collapse is insensitive to the vacancy defect density. 
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I. INTRODUCTION 

High temperature evaporation is a good method to remove impurities from carbon nan- 
otubes, typically for those samples prepared by the chemical vapor deposition method with 
Fe/MgO catalysts. This purifying technique initially found its application in the multi- 
walled carbon nanotubes (MWCNT), and has been developed to clean the double- walled 
carbon nanotubes (DWCNT) samples recently. In Ref. 1, Liu et al. apply the acid dis- 
solution and high temperature evaporation methods to purify nanotube samples. Their 
experiment shows that the residual Fe particles covered by a carbon shell can not be re- 
moved by the chemical solution, but can be cleaned through physical evaporation. A 99% 
purified sample can be achieved by evaporation treatment at 1700 °C This experiment also 
demonstrates that the samples with smaller diameter are thermally more stable than that 
of larger diameter. In Ref. |2|, Kim et al. employ the heat treatment method to investigate 
the structure transformation of bundles of single-walled carbon nanotubes (SWCNT) at dif- 
ferent temperatures. They find that the DWCNT are more stable than SWCNT and even 
comparable with the MWCNT. They also observe that outer walls of adjacent DWCNT can 
coalesce into large diameter tubes. SWCNT prefer to change into DWCNT around 2000 
°C. SHSj This is a typical phenomenon in heat treatment experiments for bundled SWCNT. 
In Ref. y, |5|, Kim et al. demonstrate that the purified bundled SWCNT evolve along very 
different pathways that depend on three different factors: initial diameter distribution, con- 
comitant tightness of the bundles, and bundle size. Samples prepared by different methods 
can result in very different final derivatives after heat treatment. These thermal evolu- 
tions are claimed to be related to various reasons, such as diffusion of carbon, C-C bond 
rearrangement, etc. 

Theoretically, the heat treatment of carbon nanotubes can be simulated by molecular 
dynamics (MD). The tight-binding MD is applied to simulate the coalesce of two ultrathin 
SWCNT after collision, and a conservation of the chiral index is found, p] In Ref. [Zl, Lopez 
et al. use the extensive MD to study the evolution of the bundles of SWCNT after heat 
treatment. They find that the bundles of SWCNT coalesce forming MWCNT through an 
interesting patching-and-tearing mechanism. 

There are two facts we can learn from above experiments and MD simulations. Firstly, 
these experiments suggest that the SWCNT are thermally unstable after 2000 °C and will 



evolve into DWCNT or MWCNT. Secondly, the underlying mechanism for the evolution of 
the SWCNT is far from clear. There are lots of possible reasons for this structure evolution. 
SWCNT samples in most experiments are in bundles and correlated by van der Waals inter- 
action. This interaction facilitates the evolution from SWCNT into DWCNT or MWCNT. 
Moreover, the residual catalyst impurities are ineluctable in SWCNT samples prepared by 
the chemical method. These catalysts may assist the evolution of SWCNT. Now a natural 
question will arise: what will happen if isolated SWCNT undergo heat treatment? How 
stable are isolated SWCNT? In a MD simulation work for the SWCNT with both open or 
closed ends, the breakdown temperature TV is found to be about 3300K, which varies in 
hundreds of Kelvin in different situations. [Sj] While in another work, the tight-binding MD 
is applied and the obtained critical temperature T^. is about 2500 K.[9] The large discrepancy 
between different existing works requires further investigations. This forms one aspect of 
the topic in this paper. 

The defects may play a significant role in the thermal stability of the SWCNT. Among 



different defects, the vacancy defect can be experimentally generated 10|, ll3|] and well 
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detected. [IJ, ll5| In Ref. Il0|, the authors apply the proton irradiation to generate vacancy 
defects in the graphite. The defect density can be quantitatively controlled by applying 
different dose of the proton. It is estimated that about 0.1% carbon atoms are displaced 
for a proton dose of 0.3 nC/fim^. The defect density can be as high as 35% for a dose 
of 250 times higher, where the carbon system is amorphous. Considering the advantage of 
these experimental techniques, it is interesting to study the thermal stability of the SWCNT 
with different vacancy defects. It has been shown that the vacancy defect can be applied 
to build a carbon nanotube semiconductor-metal intra molecular junction. [11|] The recon- 
struction and alignment of vacancies in SWCNT with different diameters are investigated 



by first-principles and tight-binding calculations. [12(] 

In this paper, the evolution of structure configuration of the SWCNT is investigated 
by MD simulation. The pure isolated SWCNT are thermally stable up to a temperature 
Tc=4655 K. The thermal vibration becomes more and more serious with the increase of 
temperature, and the cross section changes into ellipse at 2000 K. The structure deforma- 
tion of SWCNT with different vacancy defects is investigated. It is found that the critical 
temperature Tc is insensitive to the vacancy defects. This is resulting from the self-repairing 
capability of the SWCNT by heat treatment. The vacancy defect can optimize its local 



environment by exciting phonon modes localized nearby in low temperature region. At 
high temperatures, vacancy defects can migrate in the SWCNT. There are two efficient 
self-repairing processes. Firstly, different defects joint together and form a bigger hole; then 
a bottleneck junction is constructed near the defect. Secondly, local curvature is intro- 
duced to smooth the SWCNT. As a result of strong self-repairing capability, the defective 
SWCNT can find a stable configuration at high temperatures, so the critical temperature 
Tc is insensitive to the defect density. 

II. SIMULATION DETAILS 



We run MD simulation using the "General Utility Lattice Program" (GULP).[16[ The 
GULP is usually used to investigate various lattice dynamics properties of crystal systems. 
There are plenty of reliable interatomic potentials implemented in this code. This is one of 
the most important features of GULP. Another distinct feature is the use of the symmetric 
groups of lattice structures to reduce the simulation time considerably. It is worth noting 
that the MD simulation part in the GULP has been improved significantly in its latest version 
GULPS. 4. Combined with its powerful potential libraries, GULP becomes a useful platform 
to do MD simulation for the solid crystals. It is suitable for the simulation of structure 
transition and melting, as some interatomic potentials implemented can characterize the 
formation and breaking of bonds. 

We use the latest version of GULPS. 4 for linux system, freely available from its website. 
The MD simulation is run with canonical ensemble NVT at constant volume and tempera- 
ture. The constant temperature is realized by the Nose-Hoover algorithm. 
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A small 



friction parameter of 0.0005 is adopted in the Nose-Hoover to minimize the fluctuation in 
temperature. The frequency of the highest-frequency phonon mode in the SWCNT is about 
1600 cm~^, corresponding to a vibration period of about 20 fs. We use a time step of 0.001 
ps that gives 200 recording data to describe the trajectory in each vibration period of the 
highest-frequency phonon mode. Usually the total MD simulation time is 2000 ps, about 
2 X 10^ MD steps. The Newton equations are solved by leapfrog Verlet algorithm. We apply 



the Tersoff potential. 19[ This potential has a good performance in describing the formation 



and breaking of bonds. We have also tested the Brenner potential. 20j The results are qual 



itatively the same as Tersoff potential, while Brenner potential is much slower due to its 



complexity. The structure configurations in this paper are plotted by the software "Visual 



Molecular Dynamics" 2l| with trajectory data from GULP. 

Before jumping into the simulation work, it is meaningful for us to check the suitability 
of the Tersoff potential for interatomic interactions in SWCNT. Tersoff potential has found 
various applications in the description of covalent bonds. For example, this potential is 
applied to study the thermal conductance at interface between SWCNT and silicon, [22'] the 
structure of SWCNT under hydrostatic pressure, [231] and the H2 adsorption effect on the 
SWCNT. 24(] As an additional check, we use the Tersoff potential to calculate the radial 
breathing mode (RBM) in the SWCNT. On the one hand, the REM reflects the carbon- 
carbon bond strength in the SWCNT, thus is a good candidate for the check of the potential. 
On the other hand, the RBM in SWCNT has been extensively studied, so there are lots of 
existing results to compare with. RBM is widely used as a fundamental tool to clarify the 
diameter of the tubes experimentally, because it is a Raman active mode and its frequency is 



inversely proportional to t 



le tube diameter. 
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28| This mode was theoretically investigated 



by ab initio calculation, 29|] and different valence force field models. 30N32| In Fig. [H we show 
the frequency of RBM with diameter from 2.5 A to 28 A, including armchair SWCNT (n, 
n) with 2 < n < 20 and zigzag SWCNT (n, 0) with 4 < n < 36. Our results are in good 
agreement with both existing experimental data and theoretical values in whole diameter 
range. We need to stress that the Tersoff potential gives good results not only for the wider 
SWCNT but also for narrow SWCNT with diameter below 4.0 A. This shows the advantage 
of the Tersoff potential for SWCNT, because it is not an easy job to obtain correct frequency 
for RBM in narrow SWCNT, even for the ab initio calculation. Another important fact is 
that the frequency of RBM from Tersoff potential can be fitted as a function of diameter 
i^RBM = A/d where d is the tube diameter. The fitting parameter A = 2305 ± 24 cm~^A 



is very close to the first principle calculated value of 2260 cm ^A. 29j All these results 
confirm the suitability of the Tersoff potential for the description of interatomic interactions 
in SWCNT. 

III. RESULTS AND DISCUSSION 

We study the structure evolution of the zigzag SWCNT (10, 0) with the length of 26 A 
at different temperatures. There are similar phenomena observed in armchair SWCNT (10, 



10). Periodic boundary condition is applied in the axial direction. Large periods are applied 
in the two lateral directions, so that the tube is isolated equivalently in these two directions. 
Fig. [2] shows that the thermal vibration is almost invisible at 300 K, because the C-C 
covalent bonds in the SWCNT are very strong, typically in the order of 2.0 eV. The thermal 
vibration becomes obvious at 2000 K, where the cross section is not a perfect circle anymore. 
One interesting result is that the tube is far from being destroyed at 2000 K. In some 



experiments on the bundles of SWCNT, 



K 



-|5| it was found that SWCNT will always turn 



into MWCNT above 2000 K. It is possible to attribute this experimental phenomenon to van 
der Waals interaction between different SWCNT, and much more complicated environment 
in the experiments. Our simulation result indicates that an isolated SWCNT in a clean 
environment is still stable at 2000 K. However, strong thermal vibration and deformation 
can be observed in this temperature range. 

We learn from Fig. |2]that the thermal vibration becomes more and more serious with the 
increase of temperature. All carbon atoms are no longer in the same cylindrical surface at 
high temperatures. We always observe that the middle part of the tube and the two ends of 
the tube are deformed in perpendicular directions. Carbon atoms are still not evaporated 
at 4000 K, although they are in an extremely strong thermal vibration. There is no 5- or 7- 
membered rings for SWCNT observed here. Fig . |3]is the time evolution of internal pressure 



at 4000 K. The pressure is calculated by:|33l. |34| 
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where D is the degree of freedom and N is the atom number. is the interatomic interaction 
which is Tersoff potential here. This is the internal pressure and is used to balance with the 
externally applied pressure in the thermal steady state. This curve can monitor structure 
evolution of a tube. The pressure will have a rapid change if the tube changes from one 
configuration into another. Insets are some typical configurations at different simulation 
stages. 

SWCNT will be destroyed into pieces at temperatures above a critical value Tc = 4655 
K. Fig. m illustrates the collapse of SWCNT step by step. Usually the collapse of a tube 
starts from the damage of a particular region, where some atoms are evaporated. Then 
some atoms in other areas will move into this region and be evaporated again. The tube is 
finally destroyed totally. The critical temperature T^ = 4655 K is higher than two previous 



theoretical works, 



m 



9|] where the value is about 3300 K and 2500 K. One possible reason for 
this difference is the different interaction potentials which affect the value of Tc directly. In 
these two works, the total energy is either described by Brenner potential, [8| or calculated 
in the tight-binding method. [9] Another possible reason for this difference is the simulation 
time step. A large time step will accumulate numerical errors in the solution of the Newton 
equations in MD simulation. Some atoms will leave the tube solely due to the numerical 
error, leading to the collapse of the SWCNT. For example, we can artificially increase the 
time step to be 0.003 ps. It results in a value of 3600 K for the critical temperature. 
Obviously, this value is not the actual critical temperature for the tube. The boundary 
conditions in the axial direction is also very important. There are bending modes in the 
tube if both ends are fixed. These bending modes will facilitate the breakdown of the tube 



.a 



from its middle region, j 

For a longer tube with length 52 A, we find that the critical temperature is almost the 
same as the shorter one. This is because we apply the periodic boundary condition in the 
axial direction. The periodic length is large enough, so there is no much difference with the 
increase of length in the axial direction. For a wider tube (20, 0) with the same length 26 A, 
it is interesting that the T^ is significantly lower, about 4000 K. This result is qualitatively 
consistent with the experiment, |l[] where tubes in smaller diameter are thermally more stable. 



although the curvature energy in smaller tubes is larger. 35j The curvature energy is the 
internal strain energy of a SWCNT, by viewing SWCNT as rolling up a graphene sheet into 
cylindrical structure. J36|] 

Up to now, we have concentrated on pure SWCNT without any defects. From now on, 
we will switch to the discussion of the heat treatment of SWCNT with vacancy defects. 
A vacancy defect is introduced by arbitrarily removing one carbon atom from the pure 
SWCNT. Fig. [5] shows similar results as pure tubes at low temperature (300 K), where 
thermal vibration is invisible. At 2000 K, atoms around the defect will optimize their 
structure to lower the local energy introduced by the defect. These atoms have obvious 
larger vibrational amplitudes than the other atoms, which results from the localized phonon 
modes around this defect. If these localized phonon modes are excited, atoms around the 
defect have strong vibration while the other atoms almost do not move. At 4000 K, we can 
see that the defect starts to move in the tube. Fig. [6] displays how the defect shifts in a 
SWCNT. Finally, the defect will settle down in some region, and optimize the defect to an 



energy stable state. With the further increase of temperature, some atoms will be evaporated 
in the defective SWCNT but tubes are still stable. The critical temperature for SWCNT 
with only one defect is about 4635 K, which is almost the same as the corresponding pure 
SWCNT. 

To see the self-repairing process more clearly, we need to introduce more defects in the 
SWCNT. We have tried two, four, eight, sixteen, and thirty-two defects. There are totally 
240 atoms in this studied SWCNT, and the corresponding defect density for two defects is 
0.83%. It is interesting that the critical temperature Tc for these defective tubes is roughly 
the same as pure ones, around 4600 K. The fluctuation is just some tens of Kelvin. At first 
glance, the defects should result in much lower critical temperature. However, this result 
can be considered in another aspect. A stable value of Tc implies the strong self-repairing 
capability of the SWCNT. In other words, the SWCNT can repair the defects by itself and 
turns into a stable configuration. Fig. [7] shows how this self- repairing is accompanied in the 
SWCNT in case of eight defects. There are two efficient repairing methods. In the first 
process (denoted by double blue arrow), two defects move closer and closer, and finally get 
together to form a bigger hole. Then a bottleneck junction is constructed near this hole. 
It looks like a junction structure between two tubes with different diameters. In another 
repairing method (denoted by single red arrow), two defects get closer and the neighbor- 
hood around these two defects is smoothed by introducing a local curvature. Usually, the 
bottleneck junction favors to occur around a bigger defect; while the smoothing mechanism 
always happens between many smaller defects. These two mechanisms are not independent 
of each other. They can happen continuously. For example, a bottleneck junction is formed 
and immediately further smoothed. In case of large defect density, these two self-repairing 
mechanisms will cooperate with each other. We find that some defects will always stay 
along and smoothed by the heat treatment. While some defects will get together and form 
a bigger hole, which will form a bottleneck junction nearby. We want to further remark 
that the smoothing mechanism can work for self-repairing, and this repair will introduce 
a curvature surface around the defect as a byproduct. For a small hole defect, a slight 
curvature is enough to smooth the area around the defect. However, for a big hole defect, 
large curvature is in need to achieve a more stable structure by the smoothing mechanism. 
The curvature surface introduced by the smoothing mechanism is so large that it eventually 
turns into a bottleneck junction. 

8 



IV. CONCLUSION 

To conclude, MD simulation is performed to investigate the structure deformation of 
SWCNT by heat treatment. It is found that the critical temperature for the collapse of the 
pure isolated SWCNT is about 4655 K, which results from the strong covalent C-C bonding. 
With the increase of temperature, the cross section of SWCNT changes into ellipse at 2000 
K. The self-repairing processes in SWCNT with vacancy defects are studied in detail. There 
are two efficient mechanisms for the SWCNT to do the self-repairing work. The first one is 
to construct a bottleneck junction after two defects joining together. Another process is to 
introduce a local curvature to smooth the surface around the defect. Due to strong capability 
of self-repairing, SWCNT can find a stable configuration even at very high temperature. As 
a result, the critical temperature for the collapse of SWCNT is insensitive to the defect 
density. We expect that this self- repairing phenomenon can be observed in the experiment. 



13| and high 



since the generation of vacancy defects can be quantitatively controlled, 10 
vacuum chamber technique is also at a high level. 
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FIG. 1: (Color online) The frequency of the radial breathing mode in SWCNT. Results from Tersoff 
potential (this work) are compared with other theoretical or experimental values. 
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(a) T=300K 





(b) T=2000K 





(c) T=4000K 





FIG. 2: (Color online) Structure of a perfect SWCNT(10, 0) at different temperatures. Left/right 
panel are the side/top view. The colors are with respect to the position of each atom and used to 
make the figure more clear. The colors in following figures are the same. 
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FIG. 3: (Color online) Time evolution of pressure at 4000 K. The structure transition for 
SWCNT(10, 0) is illustrated as insets through some configurations of the SWCNT. 
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FIG. 4: (Color online) Time evolution of pressure at 5000 K. The structure of SWCNT(10, 0) is 
finally destroyed. Insets are some typical configurations of the SWCNT during collapse. 
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FIG. 5: (Color online) Structure of SWCNT(10, 0) with one defect at different temperatures. 
Left/right panel are the side/top view. 
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(a) t=lps 




(b) t=300ps 




(c) t=400ps 




(d) t=2000ps 




FIG. 6: (Color online) Migration of the defect (denoted by arrows) in the SWCNT at 4000 K. 
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FIG. 7: (Color online) Self-repairing in the SWCNT with eight defects at 4000 K. Two repairing 
mechanisms: (1) bottleneck formation denoted by double blue arrow; (2) smoothing mechanism 
denoted by single red arrow. 
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